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SUMMARY 

A model Is presented which can be used to Incorporate the effects of fric- 
tion and tortuosity along crack surfaces through a constitutive law applied to 
the Interface between opposing crack surfaces. The problem of a crack with a 
saw-tooth surface In an Infinite medium subjected to a far-fleld shear stress 
Is solved and the ratios of Mode-I stress Intensity to Mode-II stress Intensity 
are calculated for various coefficients of friction and material properties. 

The results show that tortuosity and friction lead to an Increase In fracture 
loads and alter the direction of crack propagation. 


INTRODUCTION 

The essential Ingredients which enter a Linear Elastic Fracture Mechanics 
(LEFM) analysis are stress Intensity factors and fracture toughness. These 
quantities can be used to predict fracture Initiation loads as well as the 
direction of crack propagation. Stress Intensity factors have been obtained 
for a large variety of loadings and specimen geometries. In most of these 
analyses cracks are modeled as discontinuities possessing smooth and friction- 
less surfaces. For Mode-I type loadings the assumption that the surfaces of a 
crack are smooth and frlctlonless Is quite satisfactory, but for Mode-II load- 
ings or mixed-mode loadings this assumption may lead to erroneous results If 
the crack surfaces are tortuous and/or offer frictional resistance to sliding. 
For such problems the effects of tortuosity and friction must be accounted for 
In order to capture the physics of the phenomena and to correctly assess the 
stress Intensity factors. This can be seen by considering the crack shown in 
figure 1(a). Let us assume that this crack was modeled as being smooth and 
frlctlonless. If a far-fleld shear stress (an apparent Mode-II loading) Is 
applied to this crack, the model would predict a zero Mode-I stress Intensity 
factor. In reality, the Mode-I stress Intensity factor Is not zero, since the 
displacement of a point along the surface of the 
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crack has both a normal and a tangential component. The tangential component 
will result In a Mode-II stress Intensity, while the normal component will 
result In a Mode-I stress Intensity. 

Cracks of the type shown In figure 1(a) have been observed In many poly- 
crystalline and aggregate materials. Some examples are metal, ceramic, con- 
crete and brittle geo- and blo-materlals such as rock and bone, respectively. 
Methods of analysis for mixed-mode fracture due to tortuosity are primitive. 
Clech et al. (ref. 1) have observed and analyzed tortuous cracks propagating 
along the Interface between cancellous bone and the cement which binds an arti- 
ficial joint to It. Their analysis considered failure under Mode-I loading, 
where the tortuosity was modeled as a no-slip Interface. The question arises 
as to what would be the most efficient way to model the problem If the struc- 
ture Is subjected to mixed-mode loading. One approach would be to model- the 
geometry of the crack surfaces explicitly In a finite-element (or boundary 
element) mesh, assign a coefficient of friction to the crack surfaces, and to 
solve the nonlinear stress analysis problem. This technique Is obviously cum- 
bersome. Concrete Is another material for which friction and tortuosity play 
a significant role during crack propagation (refs. 2 to 5). Recently, 

Ingraffea et al. (ref. 5), using the finite element method, have modeled fric- 
tional Interlock effects In concrete as distributions of normal and shear 
stresses along the Interface between opposing crack surfaces. Their model, 
however, does not account for the previously discussed coupling between crack 
sliding displacements and crack opening displacements. Riggs and Powell pres- 
ent a rough crack model for the analysis of concrete In reference 6. This 
model possesses many desirable features compared to previous attempts at rough 
crack modeling but Is complex and has a number of empirical constants that 
require specification. 

In this paper a model Is presented which can be used to characterize the 
effects of friction and tortuosity along crack surfaces n an efficient and 
realistic manner. In this model the crack surfaces are assumed to be globally 
smooth, and the roughness and friction are Incorporated through a constitutive 
law at the Interface between opposing crack surfaces. This constitutive law 
contains Implicitly the coupling between the tangential and normal components 
of the displacement of a point on the Interface (this phenomenon Is referred 
to as dllatancy), and can therefore give rise to Mode-I stress Intensity fac- 
tors even when the applied loading Is apparently Mode II. 

CONSTITUTIVE LAW 

In, this section we review the development of analytic relations between 
Increments of crack surface traction and crack sliding and opening displace- 
ments. This theory was fully developed In (ref. 7) and Is analogous to the 
theory of continuum elasto-plastlclty . The surface between two bodies coming 
Into contact Is assumed to be globally smooth as shown In figure 2(a). Atten- 
tion Is focused on a typical pair of Initially adjacent points on the contact 
surface as -shown In figure 2(b); these points are shown separated for clarity 
and the roughness that the surfaces have Is not shown and will be considered 
subsequently. 

The traction that Is supported by the crack surface at the contact point 
In question has tangential and normal components which are denoted by 
and o n , respectively. The crack sliding and opening displacements at the 
contact point pair are denoted by gt and g n , respectively, where 
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(1) 


«t = 



u 


X 


In equations (1) and (2), u x and u y are horizontal and vertical dis- 
placements (l.e., displacements tangential and normal to the global crack 
surface), respectively, while superscripts + and - denote points on the 
contact surface associated with the upper and lower bodies of material, 
respectively. 

It Is assumed that the relative displacements of the crack surface are 
addltlvely composed of reversible (elastic) and Irreversible (plastic) parts 

9 1 = 9® + gj 1 = t,n (3) 

where superscripts e and p denote the elastic and plastic parts of the 
deformation. In equation (3), the plastic part represents the permanent defor- 
mation due to sliding (and perhaps damage), which has components In the normal 
direction (due to dllatancy) as well as In the tangential direction. The elas- 
tic part represents deformations that occur Independently of sliding such as, 
for example, preslldlng deformablllty. Such recoverable deformations have been 
observed In a large variety of experiments between materials such as metals, 
ceramics, and rocks (refs. 8 to 10). The stress that the Interface supports 
Is assumed to be related to the elastic part of equation (3) by 


where the summation convention Is applied to repeated Indices and are 

the Interface stiffnesses, the determination of which will be addressed 
shortly. 


To prescribe a method for determining the plastic displacements, It Is 
necessary to assume that (1) a scalar valued slip function F = F(c»t,a n ,... 
can be defined such that F < 0 corresponds to preslldlng conditions, F = 0 
corresponds to sliding and F > 0 is not possible, and (2) there exists a 
linear relationship between Increments of traction and plastic deformation. 
Under these assumptions. It can be shown (ref. 11) that 


dg 


P . 
1 “ 



If F < 0 or dF < 0 
If F = dF = 0 


(5) 


where G Is the slip potential whose gradient gives the direction of the 
plastic deformation and X Is a nonnegative scalar that gives Its magnitude. 
Expressions for F and G can be obtained by considering the Idealized rough- 
ness shown In figure 2(c). If the friction on the active asperity surface Is 
governed by Coulomb's law, then |cq| < -yo 2 where y Is the coefficient of 
friction, a | and 03 are the tangential and normal stresses on the asperity 
surface and the convention that compressive stresses are negative Is employed. 
By stress transformation, this equation can be expressed In terms of the global 
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stresses and o n (ref. 7). This procedure leads to the following forms 

for the slip function and slip potential 

F = |<> n sin <* k + o t cos <* k | + v (o n cos « k - <* t sin <* k ) (6) 

6 = |o n sin a k + <J t COS a k | (7) 

where k = R or L depending upon which asperity surface Is active 
(fig. 2(c)). 

Combining equations (3) to (5) leads to the constitutive relation 


where 


or 


do, - dg 3 


E e P - E 

u u 


If F < 0 or dF < 0 


E e P - e 
L 1j ' U 


3£_ ^ ^ 3G_ 

3o d 1q pj 3o a 

3F E 3G_ _ „ 

3o p 3o q ' 


If F = dF = 0 


( 8 ) 

( 9 ) 


( 10 ) 


and H Is a hardening or softening parameter that can account for the damage 
of crack surface roughness. In this paper we take H equal to zero. 


The following remarks about the constitutive law should be noted: 

(1) Equation (8) Is an explicit relation between Increments of crack sur- 
face tractions and the resulting crack sliding and opening displacements. 

Unlike some other theories for dllatant crack surface problems (ref. 12) In 

which equation (8) Is postulated at the outset and the (particularly the 

E®^ and E®^) are determined either experimentally or according to some ad hoc 

rule, this theory provides an explicit and unambiguous method for determining 


(2) The theory presented In this section Is essentially a "continuum" 
theory of friction except that we refer to continuum In the sense of continuous 
area rather than continuous volume. The traction components a t and o n 
entering Into the theory are the average or macroscopic stresses that the crack 
surface supports rather than the exact stresses. Furthermore, the global crack 
surface Is smooth and the effects of roughness (e.g., dllatancy) are built Into 
the slip rule and slip potential. 

(3) The use of the decomposition In equation (3) In conjunction with 
equation (4) can be considered a relaxation of compatibility of crack surface 
displacements consisting of Impenetrability and preslldlng stick. In order to 
render the violation of compatibility Insignificant, It Is necessary to make 
the Interface stiffnesses E^j relatively large compared to the stiffness 

of the adjacent medium. In this Investigation, we take E^j = PSij, 
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where is the Kronecker delta and P is a penalty number. The penalty 

number has no physical significance and is chosen based on computational con- 
venience. Values of P that are two to four orders of magnitude greater than 
the shear modulus of the adjacent medium were used and provided solutions In 
which Incompatibilities were extremely small. Larger values of P sometimes 
gave numerical difficulties due to Ill-conditioning while smaller values of 
P gave slightly excessive Incompatibility. 

(4) It is assumed in this research that the amount of tangential sliding 
is small enough so that the asperity peak of one surface does not override that 
of the other surface. For situations in which this is not the case, the theory 
can be supplemented with an additional set of conditions as discussed In 
reference 7. 


COMPLIANCE MATRIX 


The problem which will be solved to demonstrate the significance of fric- 
tion and tortuosity is shown in figure 1. The singular Integral equations that 
govern this problem are 


2G 

ir(k +■ 1) 



f(t) dt 
t - x 


t + a, 
00 t 


-1 < X < 1 


( 11 ) 


2G f 1 e(t) dt _ 
ir(k +1) J .J t - X = °n 


-1 < x < 1 


( 12 ) 


where is the far-field shear stress, and o n are the shear 

and normal stresses applied to the crack surfaces, G is the shear modulus, 
k = 3 - 4w for plane strain, v being Poisson's ratio, and f(t) and e(t) 
are the dislocation densities, which are related to the crack sliding and crack 
opening displacements by 


- “x - J 

f(t) 

X 

dt = g t 

(13) 

-■V-J 

f 1 e(t) 

V 

dt = g n 

(14) 


Because the crack is closed at the endpoints the dislocation densities must 
satisfy the following conditions 



f(t) dt = 



dt 


= 0 


(15) 


Since and o n depend on the relative displacements between the crack 

surfaces, the singular Integral equations are nonlinear. In order to solve 
these equations efficiently, a compliance matrix will be developed that relates 
the stresses at the collocation points arising from the applied loads to the 
values of the opening (or sliding) of the crack at the integration points. The 
method relies on the numerical procedure developed by Gerasoulis (ref. 13), 
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which Is used to reduce the Integral equations to a system of algebraic equa- 
tions for unknown values of e(t) and f(t) at discrete points In the Inter- 
val [-1,1]. With the aid of the compliance matrix, It will be possible to 
obtain the crack opening and crack sliding displacements for a given stress 
distribution on the crack surfaces by solving simultaneous algebraic equations 
and hence, the need to repeatedly solve the Integral equations and Integrate 
the dislocation densities Is eliminated. The procedure for obtaining the com- 
pliance matrix Is described In detail In reference 14 and Is only reviewed 
here. Following (ref. 13), f(t) and e(t) are expressed by 


, - 1/2 , - 1/2 
f(t) -- 4>(t)(l - t^) e(t) = <l»( t) ( 1 - t') 


(16) 


and <p(t) and «p(t) are approximated as piecewise quadratic on the Inter- 
val [-1,1]. The result Is that the Integral equations are reduced to a system 
of algebraic equations through the quadrature formulas given In (ref. 13). The 
results can be written symbolically as 


C 1 E w 1 (x k )4»(t 1 ) = + ° t (x k ) k = 1 ,N - 1 


1-1 


c ! E w 1 (x k )«p(t 1 ) = * n (x k ) k = 1 ,N - 1 


£ v 4 >(t.) * 0 
1-1 1 


£ v 4 >(tJ = 0 
1-1 1 1 


(17) 

(18) 

(19) 

( 20 ) 


where c-j = 2G/ir(k * 1), N Is the number of Integration points, N - 1 Is the 
number of collocation points, and the weights w^(X|<) and v^, as well as 
the collocation points x* and Integration points t^ are given In (ref. 13) 
(the locations of the Integration and collocation points are shown schematically 
In figure 3(a)). These equations can be written as 


'l^k*! " L k 

(21) 

c l G 1k*1 * \ 

(22) 


where matrices [L] and [M] represent the nodal values of the stresses along the 
crack surfaces and the crack closure conditions. To obtain the compliance 
matrices for the crack the Inverse of matrix [G] Is obtained, and the product 

2 - 1/2 

of this matrix and (1 - t ) Is Integrated term by term to obtain a matrix 
[C] which will be called the compliance matrix for this particular geometry. 
This matrix will relate the values of crack sliding displacements (or crack 
opening displacements) at the Integration points to the values of the shear (or 
normal) stresses at the collocation points, that Is 
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(23) 


Sk( T x * w) ■ ( U J - “'x)t, - 

‘ikW ’ (“I - u y)t, ’ <9n>t, (24> 

where L Is a rectangular matrix with dimensions (N,N-1). Me wish to point 
out that 1 Because of the technique employed to determine the compliance matrix, 
the locations at which the crack opening and crack sliding displacements are 

determined (1e., the Integration points) are different from the locations of 

crack surface tractions (l.e., the collocation points). Furthermore, the 
number of displacement points Is one greater than the number of stress points. 
In order to eliminate the computational difficulties associated with differing 
numbers of collocation and Integration points, we present In the following sec- 
tion a method for rendering the compliance matrix square. 

The stress Intensity factors are proportional to the values of the func- 
tions $ and at the endpoints (ref. 13) and can be obtained by pre- 
multiplying matrices [L] and [M] by the Inverse of matrix [G] . The details of 
the calculations can be found In references 13 and 14. We note that the com- 
pliance matrix for shear loading Is the same as that for normal loading. This 

Is because the Green's functions for the two problems are the same. 


NUMERICAL SOLUTION 

Before writing equations (23) and (24) In matrix form. It Is necessary to 
make the compliance matrix square. Because the number of Integration points 
Is odd (ref. 13), there Is an Integration point at the center of the crack. 

By eliminating the equation corresponding to this point from equation (23), the 
number of displacement points Is reduced by one and the compliance matrix 
becomes square. Following this procedure the location of the displacement and 
stress points become as shown In figure 3(b). The compliance matrix Is now 
denoted by [C] and Is given by 


[C] = = 



C d * 1)J 


if 

if 


1 < ( N +• 1 )/2 
1 > (H t 1 )/2 


(25) 


where 1,j = 1,2,...,N-1. The modified forms of equations (23) and (24) can 
now be written collectively In Incremental form as 


where 




1 \ i 

1 

C I 0 
1 1 


dx + da, 
00 t 

' 


0 j c 1 
1 1 


( da n 


fd« t } 


[< d Vt, 

[ <do t»x 1 


< d9 t>t r 

X X N-lJ 

«■•*>* 1 T 

X N~1J 


(26) 

(27) 

(28) 
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with similar definitions for {dg n } and {d<* n }. d-to, Is a vector with N-l 
components, each equal to the Increment In the applied loading. 

Although the number of displacement points and stress points are now the 
same, the locations of these points are different as shown In figure 3(b). The 
constitutive law given by Eq. (8), however, relates the crack surface stresses 
at a point to the displacements at the same point. The crack displacements at 
the stress points can be expressed In terms of those at the Integration points 
by interpolation. For example, the relation between the crack sliding dis- 
placements at the stress points and the collocation points Is given by 


■ [T1 

where the coefficients of [T] are zero except for 


T 1.1 

* T 1 ,1 + 1 = °- 5 

when 

1 < 1 < (N - l)/2 

T 1.1 

= 0.75 

when 

1 = (N + l)/2 

T i,1 

+ ! = 0.25 

when 

1 = (N - l)/2 

T 1,i 

_ ! = 0.25 

when 

1 = (N * 1 )/2 

T 1,1 

= 0.75 

when 

1 = (N + l)/2 

T 1,1 

- 1 = T 1 ,1 = °- 5 

when 

(N + 1 )/2 < 1 < N 


Combining equations (8), (26), and (29) leads to 



(29) 


(30) 


(31) 


which Is a system of simultaneous nonlinear algebraic equations to be solved 
for the Incremental crack opening and sliding displacements. By taking small 
increments In far-fleld shear loading, good convergence was attained in the 
numerical work, where (E e P) was computed at each step using the displacement- 
stress configuration from the previous step and treating equation (31) as 
linear within each load increment. 


RESULTS 

The problem shown in figure 1 was solved for the various combinations of 
material properties listed In table I. Figures 4 to 6 are plots of stress 
Intensity factors as functions of the far-fleld shear stress for several values 
of the uniform asperity angle <*. It can be seen from these figures that as 
a Increases, Increases, while Ku decreases. The reason for this 
Is that frictional resistance and the kinematic constraint provided by the 
jagged surfaces resist crack sliding displacements and therefore reduce K j j . 
The dllatancy, on the other hand, produces a normal gap between the crack 
surfaces, therefore Increasing Kj. 

The curves In figures 4 to 6 have been terminated at the values of t 0 
that result In crack Initiation according to the maximum principal stress crl- 
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terlon (ref. 15), which postulates that crack growth will occur In a direction 
perpendicular to the maximum principal stress o p , when this value Is equal 
to the value of o p In an equivalent Mode-I case. The angle at which the 
crack will extend with respect to the horizontal Is given by reference 15. 



and the fracture criterion Is 

Kjc = Kj cos 3 | - 3 Kjj cos 2 | sin | (33) 

where Ki c Is the fracture toughness. 

Figures 7 to 9 are plots of the fracture stress for several materials as 
functions of a. The angle of extension e Is also shown on these figures. 

As seen from these curves, the effect of tortuosity and friction Is to Increase 
the fracture stress and to decrease the angle of extension. In order to sepa- 
rate the effects of friction from those of tortuosity, the case of steel was 
analyzed for various coefficients of friction. It can be seen from figure 7 
that friction has a significant effect on the fracture stress. When a = 40° 
and v = 0.0, the Increase In fracture stress with respect to zero tortuosity 
Is approximately 20 percent, while for y = 0.2 this Increase Is 40 percent, 
and for y = 0.5, It Is approximately 90 percent. 

DISCUSSION 

From the results obtained In this Investigation, It can be concluded that 
friction and tortuosity have a significant effect on crack tip stress fields 
when a crack Is subjected to mixed mode loading. In the example problems 
solved In this paper the Increase In fracture stress (as predicted from the 
maximum principal stress theory) with respect to zero friction and tortuosity 
ranged from 3 to 100 percent. The crack extension angles as predicted from the 
maximum principal stress theory were also affected by the Mode-I stress Inten- 
sity which arises from the dllatancy of the Interface between opposing crack 
surfaces . 

Although the use of the compliance matrix developed In section 3 renders 
the numerical problem very attractive, It Is limited to simple geometries. The 
effects of tortuosity and friction In problems with arbitrary geometries can 
be conveniently analyzed using the constitutive law in conjunction with the 
finite element method or boundary element method. 
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TABLE I. 



Kt c , 

ksl-lnl/2 

w 

V 

m 

M50 Steel 

18 


m 

m 

Concrete 

1 


mm 

Wl 

Sintered Silicon 
Carbide Ceramic 

4 


.14 

.30 
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(A) GLOBALLY SMOOTH, BUT MICROSCOPICALLY ROUGH CRACK IN AN 
INFINITE MEDIUM. 


T-oo 



(B) MATHEMATICAL IDEALIZATION. 
FIGURE 1. - PROBLEM CONFIGURATION. 



(A) GLOBALLY SMOOTH CONTACT SURFACE. 


N 

■// /)s/// 



(B) LOCAL TANGENT-NORMAL COORDINATE SYSTEM AT A CONTACT 
POINT PAIR (SURFACES ARE SHOWN SEPARATED FOR CLARITY 
AND ROUGHNESS IS NOT SHOWN). 


N 



FIGURE 2. - TWO-DIMENSIONAL CONTACT PROBLEM. 
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1 2 3 


N-2N-1 N 



1 2 N-2 N-1 

(A) LOCATIONS AS OBTAINED FROM [13]. 



1 2 N N-2 N-1 


(B) MODIFIED LOCATIONS PROPOSED IN 
SECTION 4. 

FIGURE 3. - LOCATION OF INTEGRATION AND 
COLLOCATION POINTS. 



0 50 100 150 200 


Too' 

FIGURE 4. - STRESS INTENSITY FACTORS AS FUNCTION OF 
APPLIED LOADING (STEEL). 
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Too, MPa 

FIGURE 5. - STRESS INTENSITY FACTORS AS FUNCTIONS 
OF APPLIED LOADING (CERAMIC). 
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FIGURE 7. - FRACTURE INITIATION LOADS AS 
FUNCTIONS OF ASPERITY ANaE AND COEFFI- 
CIENT OF FRICTION (STEEL). 
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MPa max 



FIGURE 8. - FRACTURE INITIATION LOAD AS 
FUNCTION OF ASPERITY ANGLE (CERAMIC). 



FIGURE 9. - FRACTURE INITIATION LOAD AS 
FUNCTION OF ASPERITY ANGLE (CONCRETE). 
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